(NASA~eR -1 20736 ) FORKING UODEL OF THE N 75 ~ 21562 ' v 

LONDON HOHENT READOUT SYSTEH Semiannual 
Progress Report^ 1 Jan« ° 30 Jmi= 1974 

(Alabama Univ ot7 Huntsville.,), 38 p HC $ 3 «75 Unclas 

CSCL 13 K G 3/34 19247 J 


WORKING MODEL 
OF THE 

LONDON MOMENT READOUT SYSTEM 


January l s 1974 - June 30 s 1974 


John Benjamin Hendricks 
and 

Gerald R» Karr 


The University of Alabama in Huntsville 
School of Graduate Studies and Research 
Huntsville 9 Alabama 35807 


SEMI-ANNUAL PROGRESS REPORT 


Contract Number NAS8-29316 


a. <.iVi 


^ <b> 


7 


Prepared for 

George C„ Marshall Space Flight Center 
National Aeronautics and Space Administration 
Marshall Space Flight Center s Alabama 35812 






\ 


March 1975 



I. PROGRESS ON STUDY PROGRAM 


This progress report covers the work done by Dr„ Karr on the theory 
of super fluid plug operation. The purpose of the work was to develop 
expressions which explain the behavior of the super fluid plug which is 
being proposed for use in liquid helium management in space for the 
relativity experiment. 

The details of the calculations and results are provided in the 
appendix. In summary, the essential results are that the mass flow rate, 
m, through the plug is given by the expression 


A 


P (T) + P, 

v n 


p S L i 
K 
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where is the total area of the plug, p is the density, S is the en- 

tropy, L is the latent heat of vaporization, and K is the total thermal 
conductivity of the plug with the liquid helium included. All properties 

are for the helium in the liquid state evaluated at the temperature on 

* 

FA 

the liquid side of the plug. The factor — - — is related to the im~ 

pedance of the pumping system on the gas side of the plug. 

Solution was also obtained for the mass flow at temperature above 
the lambda point temperature given by 




P (T) + P, 
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where p. is the viscosity,, k is the permeability of the plug and A is 
the total area of pores contained in a cross section of the plug. The 

value of A is found from measurement of the plug porosity. 

P 

Solutions were obtained for the normal and superfluid velocity 
(equation 29 and 30 of the appendix) . This was done in order to inves- 
tigate the possibility of the superfluid component reaching the critical 
velocity,, thereby violating certain assumptions in the derivation. We 
found that for the plugs of pore sizes under consideration here ? that 
the superfluid velocity is always less than the critical velocity. We 
also conclude that only for pore sizes larger than 10 micron would this 
effect become of importance. 

A brief comparison of the theory results with the experimental re- 
sults obtained by Dr 0 Urban are made. We find that the theory is very 
accurate in the region of 1.5°K but departs as much as 30% at near 2°K. 
The source of error is believed to be the approximate linear relationship 
used to represent the pumping system impedance. 
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II. PLANS FOR NEXT REPORTING PERIOD 

The next semi-annual progress report will cover work on the magne- 
tometer and low temperature research apparatus by Dr, Hendricks and work 
on magnetic shielding calculations. 




Gerald R. Karr 
Principal Investigator 


asw 
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APPENDIX 

THEORY OF SUPER FLUID PLUG OPERATION 
by 


Gerald R. Karr 


I # INTRODUCTION 

The operating characteristics of a porous plug which has liquid 
helium on one side and which is pumped on under vacuum on the other 
side is discussed. Such a device has application to the containment 
of liquid helium in a zero gravity environment. The plug has the capa- 
bility of acting as a phase separator and, to some extent, as a temperature 
control device. 

The system to be considered in the work consists of a container of 
liquid helium which is well isolated. The only means for mass flow out 
of the container is through a plug mode of porous material . The plug 
is assumed to have liquid helium on the container side while the other 
side of the plug is evacuated. In the experiment ran on earth under 
one g, the plug is situated on the bottom of the container. The physical 
system employed in the evacuation is found to be important to the plug 
operation. Three cases in particular will be considered; (1) perfect 
evacuation with zero pressure, (2) evacuation through an chocked orifice, 
and (3) evacuation through a long, small diameter pipe with heating. 

The last case is the one most likely to occur in practice. 

Two distinct ranges of temperature operation are of importance in 
describing plug flow, (1) Above T = 2.172°K and (2) Blow T = 2.172°K. 

The temperature 2.172°K is the point at which liquid helium changes from 
Hel which exists at the higher temperature range to a second liquid 
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phase Hell at the lower range of temperature. Hel is a normal liquid 

phase and the flow of this liquid is described by classical fluid flow 

relationships. In the Hell phase, however, the fluid exhibits super 

fluid properties which must be taken into consideration. 

For purposes of describing the plug operation, the two-fluid model 

of the Hell phase is used. The two-fluid model considers Hell to consist 

of a super fluid component, p , and a normal fluid component, p , where 

s n 

the total density of the fluid, p, is given by 

p = p 3 + p„ 

The ratio P n /p drops from a value of one at the X- point and asymptoticly 
approaches zero at 0°K, 

The fluid flow of Hell using the two -fluid model is described by 
the following set of Landau's equations (see for example R. J. Donnelly, 
Experimental Superfluidity , University of Chicago Press, 1967), 

o b *s - - p s P n P s - - 2 

Pa ST + P s <V 8 • v) V 8 = • — Vp + p s SVT + — V (v n -v s ) (1) 
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The first term on the right hand side represents the reaction of the 
fluid to pressure forces, the second term is associated with th® thermo 
mechanical effect, and the third term is a mutual friction dependent on 
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the motion of one component with respect to the other, Equation (2) 

contains a fourth term on the right hand side which represents the 

effect of viscosity. The flow equation for temperature above the X 

point can be obtained from equation (2) for the case of p = 0 and 

s 

p n = p. 

The flow equations are found to reduce considerably for the com- 
bination of steady state flow and small VT across the plug. These 
conditions are reasonable assumptions for plugs made with large pore 
(5-10|j) material. For small pore material (,5{J.), the VT across the 
plug is found to be large. Also, small pore materials seem to exhibit 
time dependent behavior which may require non- steady state analysis. 

In summary, the description of porous plug operation for liquid 
helium management is divided into the following items of importance. 

1. Type of evacuation 

a. perfect vacuum 

b . through orifice 

c. through pipe with heating 

2 . Temperature range 

a. above X- point 

b. below X-point 

3. Simplifying assumptions 

a. steady state 

b. small VT 

4. Velocity range of superfluid 

a. less than critical velocity 

b. greater than critical velocity 
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This report will consider the special case of small VT, steady 
state operation of a plug evacuated through an ideal orifice over the 
full temperature range. The results will be reducable to the case of 
a perfect vacuum evacuation. Future progress reports will consider the 
case of evacuation through a pipe with heating and large VT across the 
plug. 

II . SMALL VT WITH EVACUATION THROUGH CHOKED ORIFICE 

The model of the system to be described is shown in Figure 1. The 
flow through the orifice is linear in the pressure difference across the 
orifice. If the pressure downstream of the orifice is small compared 
to the pressure in the plenum, P 9 the mass flow per unit area through 
the orifice Is given by 

A/ A* = F P (3) 

P 

* 

where m is the mass flow rate, A is a flow area at the point that Mach 
Number M=l, and F is given for a choked orifice by 



where C is a orifice discharge coefficient, y is the ratio of specific 
heats, R is the gas constant of the gas and is the temperature of 
the gas in the plenum. For pressure ratios acorss the orifice of .1 


or less, the value of C is 0.85. 
s w 

In normal operation, heat Is transferred out of the container 


through the plug at a rate proportional to the conductivity across the 
plug. The heat comes from heat leaks into the container due to imperfect 
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isolation or from heat given off by experimental equipment inside the 
container. 

An analytic solution to equations (1) and (2) is obtained for the 
case of choked orifice and certain assumptions concerning the liquid 
flow in the plug. One simplification becomes apparent if equation (1) 
is added to equation (2). The result is 

P s ar + Pn aT + P s ■ V) + P n • V) 

2 - 

= - Vp + jl V v n (5) 

The right hand side of Equation (5) is the incompressible, steady- state 
Navier-Stokes equation for parallel flow. Parallel flow means that the 
flow has components in one direction only. If the left hand side of 
equation (5) can be neglected, the problem is greatly simplified. The 
first two terms of equation (5) are sero for steady-state conditions. 

The third and fourth terms are zero for the case of parallel flow and 
small temperature drop(based on continuity) „ The assumptions of parallel 
flow is justified even for flow in a porous material if one considers 
averages over regions that are large compared to the pore size. The 
continuity equation for the two~ fluid model and parallel flow is 

dp s dv g dp Sv 

v 5 s7 + P S S" + V 1 57 + Pn aT = 0 < 6 > 

Although the flow is considered incompressible, 

p + p * constant 
n n s 


( 7 ) 
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the individual densities are a strong function of temperature as 

6 

discussed earlier. For this reason, the derivatives of p and p 
may be finite if a temperature gradient exists. Since heat transfer 
is taking place, one can neglect these derivatives only if the tempera- 
ture gradient is small enough to make the resultant density gradient 
negligible. This condition can be met for plugs with large pores 
except near T^ where the density changes are large with respect to 
temperature. The large pore plugs have a high effective thermal con- 
ductivity and therefore a small gradient of temperature as will be 
discussed later. 

Under the assumption of small AT and parallel flow, equation (5) 
becomes 

2 — > 

Vp * u. V v (8) 

n 

where the velocity v^ is determined completely by the pressure drop. 

If the plug were constructed of smooth channels, such as a bundle of 
capillary tubes of radius R, the exact solution to the above for laminar 
pipe flow could be employed to find the volume flow rate given by 


4 



which shows that the flow rate is a linear function of the pressure 
drop. The plugs being considered by NASA at this time, however, are 
made of packed granular material which cannot be considered to have 
smooth passages. 

For flow through porous material, a relationship similar to 


equation (9) is available given by Darey’s law (see Muskat, Flow of 
Homogeneous Fluid®) 
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k dp 

t ► 

\i dx 


( 10 ) 


where k is called the permeability of the porous material. Values of 
k must generally be determined experimentally since reliable theory 
has not been developed for the irregularly shaped granuals contained 
in most porous substances, A good approx imat ion , however, can often 
be made by use of Koseny's equation 


1 D 3 

k «= T S- E1 -S (ID 

(1-p) s 

where p is the fractional porosity and s is the total surface area 
of the particles contained in a unit volume of the medium. Fractional 
porosity is defined as the fractional volume of voids per unit total 
volume, Thi^.vcan be obtained by weighing a sample made of granuals of 
known density. The value of s is related to the size of the granuals 
which varies considerably in most materials . 

The important result to this point is that the normal fluid velocity 
is a function of the pressure drop as given in equation (10) . For the 
experimental apparatus under consideration, the pressure on the liquid 
side of plug is the vapor pressure associated with the temperature of 
the liquid. The liquid head above the plug will also contribute to the 
pressure. The pressure drop across the plug is then given by 


- A p = P v (T) + P h - P p (12) 

where P^(T) is the vapor pressure at temperature T, is the head 
pressure and P is the plenium pressure on the downstream side of the 
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plug. The pressure P is given as a function of mass flow rate in 

P 

equation (3) . 

Considering the mass flow rate at the plug which is given by 



P v 
K s s 


+ p v 
r n n 


(13) 


where A^ is the area of pores available for flow of liquid in the plug 
and the densities as® £@e\ the liquid state. Equation (3) then becomes 


and, since 


then 


+ Pn v n 


= F 


P 
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F A* P 

P 


P 

P 



( P v 
v *s s 


+ p n V ) 
n n 


(14) 


(14a) 


Letting j£be the length of the plug, dp/dx can be written as Ap ( i . 
Equation 10 can now be written, using equation (12) and (14a) 


v 

n 




+ Pi 



(P v 


+ p 


n 



(15) 


which provides a relationship between p v and p v which will be 

3 s n n 

used later. 

In order to be consistent with the assumptions of parallel flow 
and small temperature gradient as used previously. Equation (1) should 
be written with the left hand side zero and the last term on the right 
hand side also neglected. Therefore, equation (1) becomes 



VP e p S V T 
s 
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or 

VP=pSVT (16) 

The gradient of T required in equation 16 is a function of temperature 
range (i.e., above or below T^) , the plug material, and the pore size as 
will be discussed. 

Heat is transferred across the plug in proportion to the temperature 
drop and the thermal conductivity of the plug given by 


= K V T (17) 


where A^, is the total area of the plug and where K is the thermal con- 
ductivity of the plug which has contributions from both the liquid K^, 
and the solid materials, K s given by 


K 





(18) 


( 


where A and A_ are the respective areas of the two materials. Values of 
m Li 

K as a function of temperature are available in the literature, 
m 

Donnelly discusses the effective thermal conductivity of liquid 

helium below T. , Since entropy can only be transported by the normal 
A. 

fluid component, the flow of heat then is restricted if the normal fluid 
flow is restricted. For the small openings that are encountered in porous 
materials, viscous effects strongly influence the flow of the normal fluid. 
Therefore, the smaller the flow spacing, the smaller the thermal conductivity. 
Donnelly gives as a result, for below T, 


V _ P 3 S 3 T d s 

\ ~ 12ji 


(19) 


Cl 
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where S is the entropy. Since S is a strong function of temperature, 
also becomes a strong function of temperature as evidenced by the 
approximation 


K^l.2 x 10 5 T 12 - 2 d 2 
Heat conductivity of helium at 3.3°k is only 


( 20 ) 


(3.3°k) = 6 x 10^ cal/deg/cm/sec (21) 

i 

Large pore size then results in extremely high thermal conductivities 
below which tends to be the dominant term in equation 18 . 

The heat that is transferred across the plug is liberated in 
transforming the liquid into a gas at the surface of the plug. The heat 
flow then is related to the mass flow rate and the latent heat of vapori- 
zation, L„ of the helium given by 



( 22 ) 


Equating the heat flow expression from equation (17) and (22), we obtain 


K V T 



(23) 


Substituting for VT in equation (16) using the above we get 


VP = pS 





> 


( 24 ) 


Using equation (12) and (13) 9 and (14) in equation (24) 
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+ P, 
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n n 


( 25 ) 


This equation is solved for m directly 


A 

\ 


P (T) + P. 
v h 



(26) 


Equation (26) gives the flow rate through the plug as a function of 
temperature of the liquid helium. These results are not the same as 
obtained by Selzer, Fairbank, and Everett given by 


m 


APK 

pSL 


(27) 


Equation (26) is seen to reduce to equation (27) for the special case 
of 1/F = 0 which corresponds to the case of zero pressure on the down- 
stream side of the plug. Zero pressure is, however, an unrealistic 
assumption for laboratory systems, (A factor of length, if A, is 
evidently missing from equation (27)). 


CRITICAL VELOCITY LIMITING 

Equation (26) applies only for superfluid velocities below the 

critical velocity. In order to assess the ^ladity of this solution, 

the superfluid velocity, v , must be obtained to compare with the critical 

s 

velocity, v . If the critical velocity is reached, another method of 
solution must be employed. 

In order to obtain a solution for v and v , equations (26) and 

s n 

equation (15) are solved simultaneously. Substituting for m, 


from 



equation (26) into equation (15), we get 
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(28) 


which gives the solution for v . This equation can be rewritten in 
a form similar to equation (26) 


v = 
n 


S L k 


U K 

P S L x 
K~ - 


1 

/ FA* \ 

l V 


(P V (T) 


+ V 


(29) 


the velocity of the super fluid component is obtained using equation (29) 
with equation (26) . The result is 




■r - 1- 


p p k S L A 

_JL_ 2 
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S LI 
K 


FA* 


(P V (T) + P h ) 


(30) 


The normal velocity is seen to always be a positive number while v may 

3 

become negative if 


p p k S L A 

*n r p 

K \1 Ap 


> 1 


(31) 


If the superfluid flow velocity should become greater than the critical 
velocity, equation 1 must be modified to include a viscous type term in 
the superfluid component. A solution for this case has not as yet been 
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developed since it was found that for the pore sizes and pressures under 
consideration, critical velocity was not reached. Critical velocity 
would be expected to become important, however, for larger sizes of 
pores and much higher pressure heads. 


ABOVE 

The flow rate above the transition temperature consists of only 
the normal fluid which is governed by Darey’s law as mentioned earlier, 
The flow is still effected by the orifice , however . Equation (12) for the 
pressure drop remains the same ©nd equation (14) becomes, for flow of 
vapor, 


F P = 
P 



(32) 


Using equation (10) and (12) with the above we get 

A 

P p = A*f ( P v > 


V 


JL 

\i i 


(P (T) + P, 
v h 


A p 
A*F 


v) 


(33) 


where p is the density of the liquid helium. 

Solving for p^, the flow above is given by 
Above T^ results 



( 34 ) 
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SUMMARY OF BELOW RESULTS 


Below T. 


where 


and 
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Comparison with Experiment 

The equations given above were evaluated over the temperature range 

of interest (1 • K to T.) for a number of plug parameters. The parameters 

K 

associated with plugs being tested by Dr. Urban at MSFC were employed in 
the analysis and the comparison of these analytic results with the ex- 
perimental results is presented here. The purpose of this discussion is 
to also describe the methods used to obtain some of the important para- 
meters needed in the calculation. Some of the basic properties of a few 
of the plugs to be tested by Dr 0 Urban are given in Table I, 


Table I 


Labeled 


Material 

Thickness, i 

Diameter* 

Porosity, P 

Pore 

Diameter 

Nickel 

,00630 m 

,04515 m 

.384 

10 

micron 

Stainless 

Steel 

,00630 m 

,04510 m 

,311 

10 

micron 

Copper 

,00645 m 

,04520 m 

.467 

10 

micron 

Ceramic 

,00660 m 

,04525 m 

.677 

10 

micron 

Ceramic 

,00650 m 

,04523 m 

.548 

.5 

micron 


The diameters listed are for the unmounted plug. The effective diameter 


of the plug in the holder is the same for all the plugs tested and is 


0,0437 m. 


Porosity Determination 

The values of porosity listed in Table I were arrived at by weighing 
the plug on a balance to obtain the bulk density, , of the plug given 
by 

mass 


2 „ 
TT r l 


( 35 ) 
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The porosity of the plug is then found using the density of the material, 
, in the following equation 



(36) 


The densities employed in determining the porosity were approximate values 
given in Table II, 


Table IIs Density 

of Plug Materials 

Material 

P m (gm/cm 3 ) 

Nickel 

8,902 

Stainless Steel 

7,9 

Copper 

8,96 

Ceramic s Al 2 0^ 

3, 965 


Permeability Determination 

The permeability of the porous plugs were determined in Helium gas 
flow tests ran at room temperature by Dr, Katz, The tests consisted of 
measuring the time for the plug to pass a measured volume of Helium gas 
for a measured pressure drop across the plug. The results revealed a 
linear relationship between volume of gas per unit time versus pressure 
drop over a range of pressure drops from 1/2 to 4 psi. The results of 
these measurements are given in Table III, 
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Table HI: 

Flow Rate Per Unit 

Pressure Drop 

Material 

Pore Size 

A F/A P 

Nickel 

10 Micron 

31.45 scc/sec/psi 

Ceramic 

10 Micron 

18.64 scc/sec/psi 

Ceramic 

.5 Micron 

0.448 scc/sec/psi 

Ceramic 

3 Micron 

6.23 scc/sec/psi 


The permeability of the plug is defined using Equation 10 

k = 



(37) 


Since the flow experiment measures volume flow rate rather than the 
velocity called for in Equation 37 9 multiply the right hand side by 




to obtain 



(38) 


where the term in the bracket corresponds to the measurements given in 
Table III. Taking dx to be the plug length (see Table I) s the area to 
be based on a diameter of 0.0437 m (see note with Table X), and the vis- 
cosity of the helium gas to be 

- 5 2 

jj, = 1.953 x 10 n ° sec/m 


the permeability values obtained are given in Table IV, 


Table XV; Permeability Values 


Material 


Nickel 


Ceramic 


Pore Size 


10 micron 


10 micron 


Permeability^ 


- 13 2 

3.7716 x 10 m 

- 13 2 

2.2354 x 10 m 


5.3726 x 10“ 15 m 2 


Ceramic 


.5 micron 


Effective Orifice Parameter 
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Equation 3 is a linear relationship between the flow rate and the 
back pressure developed on the gas side of the plug* This relationship 
is strictly valid only for a sharp edged orifice. The pumping line em- 
ployed in the experimental apparatus to be considered here has a minimum 
diameter of 1/4 inch with many turns and a length of over six feet before 
expanding to a larger diameter. The pumping line reaches room temperature 
soon after leaving the dewar. 

The restriction to the gas flow caused by the small diameter, long 
tube with high heating is represented in this analysis as a single sharp 
edged orifice. Modeling the pumping line restriction as a sharp edged 
orifice is not particularly accurate but serves to represent an important 
part of the plug system in a simple manner. Future work will incorporate 
a more accurate representation of the pumping line restriction which calls 
for an equation of the form, 

m = G P p 2 (39) 

instead of the linear relationship given in Equation 3. The F factor in 
Equation 3 is a function of temperature on the gas side of the plug as 
given in Equation 4. This temperature dependence will be ignored in the 
present treatment since the error committed is less than the error of using 
Equation 3 rather than Equation 39. 

The determination of F and A cannot be made separately but an 
effective orifice parameter can be found by writing Equation 3 as 



( 40 ) 
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In Equation 40, the values of m, A,^ and P^ are measurable quantities 
allowing the determination of the effective orifice parameter, placed in 
brackets in Equation 40 a In one of the first plug experiments at MSFC, 
the values of m, A,^ and were monitored in order to obtain a value 
of FA / A,^o From these measurements the following value was obtained 



1 O 0412 x 10 


JSg 

m „ sec 


(41) 


which is representative of the pumping line restriction at temperature 
near 1 D 6 °K where the pumping rate was found to be near 500 scc/min 0 

Liquid Helium Properties 

The following liquid helium properties are required for the com- 
putation of plug flow rates „ 


S 

L 

Pn 

Ps 

P 

M 


Vapor pressure 
Entropy 

Lantent heat of vaporization 
Density of normal fluid component 
Density of super fluid component 
Total density of fluid = P s + P n 
Viscosity of the normal fluid component 


The vapor pressure is obtained from a curve fit relationship provided by 

2 

NBS which gives the vapor pressure in units of n/m 


P v = 1 0 33322 x 10 5 EXP a l T + a 2 + a 3 /T + a^/T 2 


+ a l4 /T 


12 


(42) 
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where 


1 

= 

-49,510540356 

2 

= 

651,92364170 

3 

= 

-3707,5430856 

4 

= 

12,880,673491 

5 

= 

-30,048,545554 

l 6 

as 

49,532,267436 

l 7 

= 

-59,337,558548 

l 8 

= 

52,311,296025 

l 9 

- 

-33,950,233134 

l 10 

= 

16,028,674003 

L 11 

= 

-5354,1038967 

l 12 

- 

1199,031906 

l 13 

= 

-161,46362959 

l 14 


9,8811553386 




The other liquid helium properties were obtained from 12 point orthogonal 
polynomial least squares curve fits of values given in Donnelly, Experimental 
Superfluidity ^ 1967 „ The values used are given in Table V. The density 
of the super fluid, p and the total density, p, are obtained from the 
values of p^ and the ratio of p^ to p given in the table. Therefore, 


P s = P n ' (p n /p) ' p 


(43) 



Table Vi Liquid Helium Properties 


T 

Temperature 

°K 

S 

Entropy 
j / gm°K 

L 

Latent Heat 
of Vaporization 

cal/mole 

^n 

Normal Fluid 
Density 

gm/ cnf* 

p n ; p 

Ratio of Normal to 
Total Fluid Density 

Id 

Viscosity of 
Normal Fluid 
Micro-poise 

1.1 

O 0304 

19.75 

.00212 

.0146 

26.8 

1.2 

O 0523 

20.30 

.00405 

.0279 

17.5 

1.3 

.0853 

20.80 

.00686 

.0473 

15.2 

1.4 

.132 

21.20 

,01100 

.0743 

14.1 

1.5 

„ 197 

21,48 

.01650 

.1140 

13.5 

1.6 

a 284 

21,78 

.02465 

.1700 

13.2 

1.7 

.398 

22,05 

.03460 

.2380 

13.0 

1.8 

.535 

22.20 

,04720 

.3250 

13.0 

1.9 

.709 

22.26 

.06240 

.4290 

13.4 

2 o 0 

.929 

22.18 

,08070 

.5540 

14.9 

2.1 

1.215 

22.02 

,10550 

,7230 

18.6 

2,2 

1,600 

21.87 

.14600 

1.000 

28,0 


ro 

Ul 



and 
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P = 


P„ + P, 


(44) 


Thermal Conductivity 

Values of the thermal conductivity of the plug materials were ob- 
tained from National Bureau of Standards, NBS Monograph 131, Thermal 
Conductivity of Solids at Room Temperature and Below, 1973 „ At the 
temperatures under consideration, the material thermal conductivities 
are found to vary in a linear manner on a log- log plot of conductivity 
and temperature „ The thermal conductivity is then written as 


K m = EXP (A + B log T) 


(45) 


where A and B are the parameters of the straight line curve fit. 

Thermal conductivity is found to be a strong function of the im- 
purities or the alloy of the material under condideration. Table VI gives 
a sampling of the thermal conductivities of various materials at 2°K 0 
Figure 1 shows the variation of thermal conductivity with respect to 
temperature for a number of materials. 



Figure 1: 


Thermal conductivity variations of mernls and 

alloys f — ), electrical insulators ( ) 

and a semiconductor (o — o_ 
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Table Vis 


Thermal Conductivity of Materials at 2 K 


Material 


K m s * ° K ) 


Nickel -high purity RRR - 670 

-Matthey 37043, 99.995 pure 
- 0 16 atomic percent Cu 

Silver -better than 99.99% pure 


Gold -99,999% pure with traces of Cu, 

Ag, Si, and Pb . RRR = 536 

Copper -.003% Ag, Ni, and Pb ^ 

-Commercial wire RRR * 85 


Stainless 

Steel -Extrapolation of Mean of number 
of stainless steels 


Alumina 

A1 2 0 3 


-Sintered to within 5% of crystal 
density, estrapolated value 
-Doped with .1 (weight) 


Fused Quartz 


-Extrapolated from 100 K 


Aluminum -99,998% pure 

-Alloys (extrapolated) 

Titanium -99,99% pure 

-98 % pure 


550.0 
47,5 
20.0 

370.0 


125.0 

118.0 
260,0 


0,10 


.3 

13.5 


.08 

480.0 

2.0 


1.0 

.4 


For comparison, Table VII gives the effective thermal conductivity of 
liquid helium at 2°K calculated using Equation 19, 


Table VII: Thermal Conductivity of Liquid Helium at 2°K 

(Equation 19) 


Pore Size (microns) of 
flow passage 


K 

L 


W 


m 


K 


0.1 

0.5 

1.0 

5.0 

10.0 


20.46 

511.54 

2,046,2 

51,154.0 

204,620.0 
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The thermal conductivity of the liquid is seen to be generally higher 
than the material values given in Table VI except for pore sizes less 
than ,5 micron in diameter 0 

The thermal conductivities of the porous plug materials are ex- 
pected to be less than the bulk material because of the granular structure 
of the plug. Since thermal conductivity measurements were not made of 
the porous materials used in the plug tests, the bulk properties are em- 
ployed in the analysis with the realization that these conductivities 
are likely higher than that of the porous material 0 The error made in 
using the higher conductivity values is not important for large pore 
plugs since the liquid helium thermal conductivity dominates the total 
conductivity, For cases in which the conductivity of the plug material 
is important (for pore sizes less than 0,5 micron), better values of 
thermal conductivity are needed in order to accurately predict the plug 
operation. While the relationships used here are expected to be in error 
in absolute value, the variation with respect to temperature is repre- 
sented, The characteristics of small-pore-size plugs revealed in ex- 
periments may yield improved values of thermal conductivity when such 
results are compared with those predicted. 

In view of the lack of porous material conductivity information as 
discussed above, curve fit parameters A and B were derived to represent 
a number of different materials. These values are given in Table VIII, 
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Table VIII; Curve Fit Parameters 

Material Represented 

Nickel - high purity 
Aluminum Alloy 
Copper 

Alumina (Al^ 0^) 

Beryllum Oxide (B e 0) 

Stainless Steel 
Quartz 

Values of A and B listed here, when 
conductivity in units of W/em°K. 


* 

for Thermal Conductivity 


A B K (2°K) 




- . 17712 

.78873 

144.7135 

-4.8655 

1.0269 

1.5706 

.2332 

1.008 

253.93 

-3.141 

1.348 

11. 

-9.388 

2.817 

,059 

- .8011 

1.239 

105,94 

-8.087 

.8472 

.055 

used in Equation 

45, will give 

thermal 


Area of PoreSt, Ap. and Area of Material,, Am 

The effective area available for flow of liquid helium is obtained 
in the following derivation,. Let the flow passages for the liquid helium 
be approximated by straight tubes of constant diameter 0 If each tube 
were to have the diameter of the pore size of the plug s then each tube 
would have an empty volume of 


Volume of one flow tube = 



(46) 


If it is assumed that the flow tubes comprise all of the empty space of 
the plug s the total volume of voids, V 9 is 


W/m °K 
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(47) 


where N is the number of flow tubes . The value of N for a particular 
plug is found from the porosity measurements as follows. Since porosity 
is defined as 



(48) 


where is the bulk volume of the plug. The value of N is then found 
from Equations 47 and 48 



(49) 


The total area of pores is found to be 


A 

P 



P V 
^b 

Z 


P 


Z 


(50) 


The following expressions are therefore obtained 

A p = P A t (51) 

and 

A„ = K,, - Ap = (1-P) A^, (52) 

These estimates of and A^, will be employed in the relationships de- 
scribing the liquid helium flow in porous materials. 
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Pore Size Estimation 

The size of the pores in the porous plugs tested are not uniform as 
evidenced by electron microscope pictures made of a number of the plugs. 
The labeled values of pore size are very approximate and is meant as only 
an estimate of the pore size. Since pore size is required in the eval- 
uation of the equations presented here, it was decided to relate the 
pore size to the measured value of permeability rather than use the 
labeled value. 

The permeability measurements given in Table III were least squares 
fit with respect to the labeled pore size. A linear curve fit was em- 
ployed with the result 


d = ( A F / A P) / 2.5 (53) 

where d is the pore diameter (assumed to have a circular cross section) 
and A-F/A P corresponds to the room temperature helium gas flow given 
in Table III in units of standard cubic centimeters per second per psi 
of pressure drop (SCC/SEC/psi) . The value of d determined using 
Equation 46 will be given in units of microns. 

Results 

The above relationships were evaluated on the UAH Univac 1108 com- 
puter. The computations were performed using the international system 
of units. The flow rates were then converted to units of Standard Cubic 
Centimeters per Minute (SCC/M) for ready comparison with the experimental 
measurements which were made in those units. Table IX gives the results 
for the first plug tested which was nickel having the properties listed 
with the table. Table X gives results of the 10 micron ceramic and Table 
XI gives results for ,5 micron ceramic. 
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Table IX; Results for Nickel Plug (10 Micron) 

Material = Nickel Porosity = .38411 

Pore Size = 12.58 Thickness = .0063 m 

Thermal Conductivity at 2°K = 144 W/m°K 


Permeability 

- 13 2 

= 3.77 x 10 u in 

Temperature 

°K 

Flow Rate 

scc/m 

1.1 

20.3 

1.2 

43.6 

1.3 

84.4 

1.4 

150.7 

1.5 

251.7 

1.6 

398.1 

1.7 

601.0 

1.8 

872.3 

1.9 

1,223 . 3 

2.0 

1,663.8 

2.1 

2, 200.0 

2.16 

2,568.0 


m/ sec m / sec 


.538 

X 

h-> 

o 

1 

Ul 

.653 

X 

io “ 6 

.116 

X 

10" 4 

.127 

X 

10- 5 

.152 

X 

o 

I— t 

.237 

X 

10" 5 

.176 

X 

10- 4 

.432 

X 

10" 5 

.191 

X 

10" 4 

.761 

X 

10" 5 

.201 

X 

-4 

10 

.130 

X 

10" 4 

.208 

X 

10“ 4 

.216 

X 

10" 4 

.213 

X 

o 

1 

.356 

X 

- 4 

10 

.214 

X 

-4 

10 

.599 

X 

10“ 4 

.210 

X 

10’ 4 

.106 

X 

10" 3 

.201 

X 

10‘ 4 

.229 

X 

10" 3 

.192 

X 

10" 4 

.563 

X 

10" 3 


i 
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Table X; Results for Ceramic Plug (10 micron) 

i 

Material = Ceramic i Porosity = .677 

Pore Size = 7.456 micron 1 Thickness = .0066 m 

Thermal Conductivity at 2°K = 11 W/m °K 

Permeability = 2.2354 x 10 ^ 3 m^ 

Temperature Flow Rate V g 

°K SCC/M m / sec m/ sec 


1.1 

20.0 

.113 

X 

1 

o 

r*~t 

.241 

X 

io“ 6 

1.2 

43.4 

,142 

X 

-4 

10 

.496 

X 

10" 6 

1.3 

84,2 

.158 

X 

10“ 4 

.986 

X 

10“ 6 

1.4 

150.5 

.174 

X 

o 

.-t 

.185 

X 

io“ 5 

1.5 

251.6 

.186 

X 

10" 4 

.332 

X 

10“ 5 

1.6 

397.9 

.194 

X 

— 4 

10 

,573 

X 

io“ 5 

1.7 

600.8 

.200 

X 

10“ 4 

.966 

X 

10“ 5 

1.8 

872.1 

.204 

X 

!0- 4 

s 162 

X 

-4 
10 H 

1.9 

l p 223 .0 

.205 

X 

10“ 4 

.277 

X 

-4 

10 

2.0 

l s 663 , 6 

.202 

X 

ID' 4 

.500 

X 

-4 

10 

2,1 

2 S 199 . 7 

.192 

X 

-4 

10 

.109 

X 

10“ 3 

2.16 

2p 567 . 7 

.184 

X 

i 

o 

r— f 

.270 

X 

10" 3 




1,1 

16 ,0 

.260 

X 

io' 5 

.366 

X 

10‘ 6 

1,2 

30,7 

, .117 

X 

io" 4 

.455 

X 

io" 6 

1.3 

52,2 

| ,336 

X 

io -3 

-.304 

X 

io" 6 

1.4 

83.9 

.747 

X 

b-> 

o 

i 

■P* 

-.377 

X 

!0- 5 

1.5 

130,9 

,141 

X 

io" 3 

-.146 

X 

io’ 4 

lo 6 

. 204.5 

,231 

X 

io' 3 

-.409 

X 

io" 4 

1.7 

321.7 

.338 

X 

io -3 

-.954 

X 

io" 4 

1,8 

502.7 

.448 

X 

io- 3 

-.197 

X 

io" 3 

1,9 

761,3 

.543 

X 

io -3 

-.374 

X 

io" 3 

2.0 

1,096 . 1 

,600 

X 

io- 3 

-.685 

X 

io -3 

2,1 

1,484.9 

,606 

X 

io" 3 

-.145 

X 

io- 2 

2,16 

1,717.4 

.580 

X 

io“ 3 

-.351 

X 

io" 2 
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Discussion of Results 

While a large number of tables could be generated, the three tables 
of results provided here are sufficient to illustrate many of the basic 
characteristics of all the results obtained using the relationships 
provided , 

a. Comparison with experiment 

Table XXI gives samples of the results obtained in experiments 
conducted at MSFC and Stanford. 

Table XII s Experimental Results 


Temperature 

MSFC 

Nickel Plug 
7-10 micron 
SCC/M 

Stanford 
Rolled Aluminum 
Foil Plug 
SCC/M 

1.6 

3 70, 

312, 

1.7 

540. 

564. 

1.8 

760. 

803. 

1.9 

1200. 

624. 

2.0 

1900, 

442. 

2.1 

3000, 

336. 

2.16 

3800. 

336. 


The Stanford results were obtained from "Research at Stanford on the 
Containment of Liquid Helium in Space by a Porous Plug and a Long Hold- 
Time Dewar for the Gyro Relativity Experiment" by Lipa s Everitt, and 
Fairbank „ The fact that the MSFC plug and the Stanford plug had 
different dimensions was not taken into consideration in constructing 


Table XII . 



The experimental results given in Table XIX show similar total flow 


rates in the range of 1,6 °K to 1,8 °K for both the MSFC and Stanford 
plugs. However, at 1 0 8 °K and higher, the Stanford plug shows a reduction 
in flow rate with increasing temperature to the lambda point. The MSFC 
plug on the other hand shows a monatomic ally increasing flow rate with 
increasing temperature which is in agreement with that predicted by the 
equations developed here as given in Tables IX, X, and XI, 

The theory developed here is found to predict the behavior of the 
MSFC large pore Nickel plug with respect to temperature. While the 
shape of the floes rate versus temperature curve is accurately predicted, 
the values predicted are about 307, , less than that measured at temperatures 
near the lambda point. The agreement at lower temperatures is progressively 
better until the best agreement is reached at 1,6 °K, (7,5% error). 

The reason the theory has such large error at the higher temperatures 


was traced to the value of the pumping factor 


FA* 


used in the program. 


This factor was determined for flow measurements made near 1.6 K. If 

I 

the pumping line restriction were an ideal choked orifice as modeled, 
this factor would be constant, Hcwever, since the restriction is more 
complex, the factor is not constant and is found to increase as the 


pressure on the gas side of the plug increases. Future work will modify 
the model of the pumping line restriction to yield improved comparison 
with experiment. 
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